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Abstract 

The Radiotherapy treatment planning optimization process based on a quasi-Newton algorithm with 
an object function containing dose-volume constraints is not guaranteed to converge when the dose 
value in the dose- volume constraint is a critical value of the dose distribution. This is caused by finite 
differentiability of the dose-volume histogram at such values. A closer look near such values reveals 
that convergence is most likely not at stake, but it might be slowed down. 

1 Introduction 

A recent development in Radiotherapy treatment planning is to use an optimization process to obtain a 
treatment plan. To do so one needs to formalize the clinical objectives in a so called object function whose 
optimum should correspond to an optimal treatment plan. Many current treatment planning systems rely 
on a quasi-Newton algorithm to perform the optimization. In order to guarantee that this process indeed 
converges to an optimum the object function should at least be twice continuously differentiable. This 
last property is the main subject of this study. As it turns out the object functions used in Radiotherapy 
do not always have this property so that convergence of the optimization process can not directly be 
guaranteed in those cases. 

In section [2] we present some background for the use of optimization in Radiotherapy treatment planning. 
The object functions currently used are built from so called dose- volume and EUD constraints. By dose- 
volume constraints, the dose-volume histogram enters the object function. It is the finite differentiability 
of the former that causes the insufficient differentiability of the latter. In section EH this relationship is 
explained in more detail. A example illustrates how finite differentiability of the object function influences 
the optimization process. The practical consequences thereof are the subject of section [4] Where section 
[3]is kept informal, section [5] provides formal statements and proofs of the differentiability results. This 
last section does not contain new arguments but serves as a support for the exposition. 

2 Background for optimization in Radiotherapy 
2.1 Optimization in Radiotherapy 

In Radiotherapy treatment planning a careful choice of values for treatment parameters is necessary to 
obtain a dose distribution in the patient that meets criteria for eradicating a malignant tumour. At the 
same time it should meet criteria so that healthy tissues and organs are not critically damaged. Due 
to the fact that dose is deposited in every irradiated part of the patient, these criteria are conflicting in 
almost every clinical situation. Experienced planners are usually very good at intuitively balancing such 
conflicting goals, still there is a need for more objective and quantifiable methods to devise and evaluate 
treatment plans. This has led to the construction of object functions, depending on treatment parameters, 
which formalize the objectives of dose on tumours and healthy tissues, see [U [El |9]. By optimizing this 
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object function with respect to treatment parameters, one hopes to get an optimal treatment plan. 
Another, very different, reason to use optimization in treatment planning is the emergence of IMRT 
where the sheer number of treatment parameters simply precludes manual planning. 

Even if we know that a function has a unique optimum, there is no simple recipe to find it. Therefore 
we have to resort to numerical, iterative methods to solve optimization problems, which, given a starting 
point, generate a sequence converging to an optimum. The latter can be guaranteed only if certain 
conditions on the object function are met. It seems, however, that if these conditions are considered at 
all, it is taken for granted that they are fulfilled. 

Many treatment planning systems use a quasi-Newton method, like the LBFGS algorithm, to find an 
optimum of the object function, see for example [5] or [7] as a general reference. One of the conditions 
guaranteeing that this algorithm yields a sequence locally converging to an optimum is that the object 
function is twice continuously differentiable. The object functions currently used in treatment planning 
systems are limited to so called EUD constraints and dose-volume constraints, or combinations of both, 
see for example [8] , but also see section 12.21 remark [TJ Our results show that the former do satisfy 
the differentiability condition. The latter however, do not, that is not for all values of the treatment 
parameters. 

2.2 Object functions in Radiotherapy 

The goal of Radiotherapy treatment planning is to devise a treatment plan achieving a high tumour 
control probability (TCP) combined with a low normal tissue complication probability (NTCP). A priori 
both TCP and NTCP are unrelated to treatment parameters. Even relating them to criteria on the 
dose distribution in the patient is a non-trivial task. To simplify matters one does not consider the 3- 
dimensional dose distribution, but a derived quantity namely the dose- volume histogram for each relevant 
region. The latter are tumours and healthy organs for which the dose- volume histogram tells which part 
of the volume receives at least a certain dose, see equation ([T]) below. Data from previously treated 
patients are used to relate TCP and NTCP to dose- volume histograms, see [3]. At this point further 
simplifications are necessary, because TCP and NTCP can only be related to certain points of dose- 
volume histograms or even only to averages. Generally for tumours this leads to dose- volume constraints, 
that is a minimum and a maximum dose for volumes t> m in and w max respectively or in a more concise 
notation (d m in, ^min) and (d max ,v max ). In case of healthy tissues the situation is much more complicated, 
in practice one uses a mixture of dose- volume constraints and EUD constraints, see [8]. 

Criteria for the dose distribution or the derived dose- volume histograms are related to treatment param- 
eters, thus the treatment planning objectives can be expressed with an object function of these treatment 
parameters. Now the object function used in treatment planning is a weighted sum of 'local' object func- 
tions, one for each dose- volume constraint or EUD constraint. Let us consider a tumour with a minimum 
and a maximum dose- volume constraint. Furthermore V c is the dose- volume histogram of the tumour 
for a given dose distribution and treatment parameters a, that is V a (h) is the (relative) volume of the 
tumour, region R, receiving at least dose h 



In this expression vol(A) is the volume of the region A and j a is the dose distribution. One constructs 
what we will call local object functions. For example the local object function for a minimum does- volume 
constraint for a specific volume is a function i^ m i n (c) which increases the more the actual relative volume 
at d m i n is below v m i n and decreases or even becomes zero if the actual relative volume at e? m i n is above 
Vmin- Similarly a maximum dose object function i 7, ma x(c r ) is constructed such that it increases when the 
relative volume at <i max is above f max and decreases or becomes zero otherwise. Thus the object function 
for the tumour alone would be the weighted sum of two local object functions, usually of the following 
form 

F(a) = Wi F min (a) + w 2 F max (a) = wi G(v mlu - V^drnin)) + w 2 G(V a {d max ) - v max ), (2) 
where w\ > and W2 > are weights and 



V a {h) 



vo\{x £R \ fa(x)>h}. 



(1) 




(3) 
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Although the object function F is a function of the treatment parameters a it does not depend on them 
directly, but only via the dose-volume histograms. The latter depend on the dose distribution which 
directly depends on the treatment parameters. 

The EUD constraints are in a similar way expressed in the equivalent uniform dose. A quantity used to 
measure the biological effect of radiation especially for healthy tissues, see [6]. We denote the equivalent 
uniform dose by E a , here it is a function of treatment parameters a 



1. The word constraint in this context is somewhat misleading because a dose- volume nor an EUD 
constraint is used as a constraint in the sense of constrained optimization. If they were, many if 
not most optimization problems in treatment planning would have no solution. It is a common 
clinical situation that a tumour is located adjacent to a critical organ. For several reasons the 
region, called PTV, on which the treatment dose is prescribed is larger than the tumour. See [5] for 
a systematic description of XTV's, where X = G, C, I or P. Therefore the critical organ is possibly 
partly located in the PTV. Suppose for example that one third of the volume of the critical organ 
is located inside the PTV. Furthermore suppose that the minimal dose-volume constraint on the 
PTV is 60 Gy (Gy is unit of dose) and the maximum average dose constraint on the organ is 20 Gy. 
Since the dose distribution varies continuously in space and its gradient is bounded, this constrained 
optimization problem will have no solution. Due to the anatomy a minimum dose of 60 Gy on the 
PTV implies an average dose well above 20 Gy on the organ and vice versa. A weighted sum of 
local object functions is a way to find a clinically acceptable compromise for this situation. To 
avoid the term constraint some authors use the term soft constraint, but in my opinion this term is 
superfluous because we already have the notion of object function. Nevertheless we adhere to the 
word constraint, in order to stay in line with the existing literature. 

2. The equivalent uniform dose, as defined in equation (01, can be regarded as a generalized average, 
in fact for a = 1 it is the standard average of the dose distribution over the region R. Formally it 
is identical to the L p -norm of the dose distribution f a restricted to a region R of the patient. But 
p is called a and it has a biological interpretation. One associates a value of a to each tissue type, 
but there is no agreement yet in the literature about these values. 

3. The dose-volume histogram refers to the graph of a function that is usually given the same name. 
This function, defined for each a in equation ([T]), can be regarded as a volume function in the 
sense of [4]. Let us fix a for the moment. The dose- volume histogram of a region R is a function 
h i ► V a (ti) whose value is the volume of R receiving at least dose h. In other words it gives the 
volume of a region enclosed by two level sets of the dose distribution f a , namely those x £ R where 
f a {x) = h and the x G R where f a (x) — maxn(f a ) (this last level set will be a point). 

4. Any quantity that influences the dose distribution may be considered as a treatment parameter. 
Usually treatment parameters are beam intensities and beam angles, that is parameters of the 
treatment unit. 

5. The object function F(a) in equation ^ also depends on d m i n , d max , v m - ln , u max , w\ and W2- Such 
parameters might be called planning parameters. Here we will keep them fixed. 



3 Differentiability of object functions in Radiotherapy 

The numerical optimization process based on a quasi-Newton method converges locally to an optimum 
if the object function is at least twice continuously differentiable. In the sequel we will simply say 
diffcrcntiable. In current Radiotherapy practice the object function is a function of treatment parameters, 
but only indirectly via dose- volume constraints and or EUD constraints. Here we will consider their 
differentiability properties. 




(4) 



Remarks 
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Let us start with object functions only containing EUD constraints. Then F(a) is a weighted sum of 
local object functions involving only 

G(d m i n — E a (a)) and G(E a (a) — 

like in equations @ and ©. In this case the object function F is differentiable. Here we have to assume 
that the dose distribution is differentiable. 

However, when the object function contains dose- volume constraints it is not always differentiable. This is 
due to the finite differentiability of the dose-volume histogram at critical values of the dose distribution, 
even when the latter is differentiable. Since it is not immediately obvious how the differentiability of 
the dose-volume histogram as a function of dose affects the differentiability of the object function as 
a function of treatment parameters, we present the following explanation. Essentially, the treatment 
parameters enter the argument of the volume function because in the presence of parameters, the critical 
values are parameter dependent. 

When treatment parameters vary, the dose distribution f a varies and in particular its critical values. 
Thus a critical value may pass through h. Suppose x a — is a critical point for a — a® (a critical point 
can always be translated to 0). It can be shown, see section El that for values of a near Co and x near 
0, f a can be transformed to a new function f a (x) = g{x) + f a {x a ). Where g is a a independent local 
standard of / CT with critical point at zero and critical value zero and f a {x a ) is the critical value of f a for 
a near a®. 

Using the previous we will show in section \E\ that V a (h) = lgiob(c, h) + Vi oc (a, h) can be split into a global 
and a local part. The global part depends differentiably on <r. Let us consider the local part in more 
detail. On a neig hbourhood O C R 3 of we have 



The function V\ oc (a,h) is proportional to the volume in the expression above and the proportionality 
factor is a differentiable function of a. The last line in the formula above shows the volume function 
V g for a region enclosed by level sets of g. According to [I] this function is not twice continuously 
differentiable at critical values of g. Thus Vi oc (a, h) is not twice continuously differentiable at values of a 
for which 



This, in turn, means that the object function F is not twice continuously differentiable at values of a 
such that h is a critical value of f a . 

Remarks 

1. The function G in equation ((3]) is not differentiable at 0. This alone would already imply that the 
object function in ^ is not differentiable. However we will assume that we are in a situation like 
in remark [T] where we are looking for a compromise. This means we stay away from points where 
G is not differentiable. 

2. Note that equation ([SJ defines a hyper surface A in parameter space M. m . It will be to— 1-dimensional 
in most of its points, therefore we need to vary only one parameter to pass through A. 

3. Throughout we assume that the dose distribution is a differentiable function. If it were less differ- 
entiable, the volume function would also be less differentiable for both critical and regular values. 

3.1 Examples 



V loc {a, h) cx vol{a; e O \ f a (x) > h} 



= vo\{x E O | g(x) + f a {x a ) > h} 
= vo\{x g O | g(x) >h-f a {x a )} 
= V g (h-f (T (x (T )). 



(5) 



The aim of this section is to provide two examples. One shows a differentiable object function, the 
other shows an object function that fails to be twice continuously differentiable at a certain value of 
the treatment parameters, here we take a = 0. In the first case the Newton iteration converges to 
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an optimum, wheras in the second case the Newton iteration does not converge if one of the iterates 
becomes 0. To avoid lengthy calculations obscuring the points we wish to illustrate we will make several 
simplifications. 




x x 
Figure 1: On the left the graph of f ai (x, 0, 0) and on the right the graph of f a2 (x, 0, 0). d + , d max and d c 
are explained in the text. 



Let f a be a dose distribution on the region R. The object function consists of a single maximum dose- 
volume constraint (<i max , %ax) so F(a) — G(V a (d max ), v ma x)- The difference between the to examples will 
be the dose distribution. Suppose for a in a sufficiently large neighbourhood of <J\ the dose distribution 
is a function with a single maximum and for a similar neighbourhood of a% it is a function with two 
maxima for example: 

f ai (x,y,z) = — — 1 — r and f a2 (x,y,z) 



l + x 2 + y 2 + z 2 3V ' l + x 2 + y 2 + z 2 2 + {x - A) 2 + y 2 + z 2 ' 

Here o\ and 02 are different parameter values. In figure [T] we draw the functions f ai (x, 0, 0) and 
fa 2 (x, 0, 0). For both dose distributions we call d+ the value at the global maximum and d c is the 
value of the other local maximum of / CT2 . Both d + and d c depend on the treatment parameters cr. Let 
us assume that always d c < d max < d+. Furthermore we assume that Va(d max ) < V a (d c ). 

For cr near <j\ the dose distribution has only one critical value namely d + so for d < d + , V a {d) is 
a differentiable function of a which we call U. In the second case, where a is near 02, we will be 
particularly interested in d = d c . For d near d c but d ^ d c , V a (d) is again a differentiable function. At 
d — d c , V a (d) = U(a) + V g {d c — d) can be split in a local and a global part, as in the previous section, 
where the global part is differentiable. We assume for simplicity that the global part is again U. The 
local part V g is essentially a function of the following form 



V g (h) 



ah? , h > 
0, h<0, 



where g is a local standard form of the dose distribution, a is a constant. Since the volume function V g is 
not differentiable at 0, the object function (that we will construct shortlye hereafter) is not differentiable at 
values of cr for which d = d c . The latter defines a hyper surface A in parameter space of codimension one, 
which means that we need to vary only one parameter to cross A. Therefore we choose new parameters 
and the first one we take to be d — d c . Only this parameter will be relevant for the example, we call 
it again cr, but now a 6 M. Thus for the second dose distribution we finally have the following volume 
function: 

v /„ _ / U(a) +a(-a)i, a<0 

Va{ ) ~ \ U(a), a > 0. [b) 



Now we construct the object function. First we scale such that w max = 1. Then the object function for 
the first dose distribution is 

F 1 {a) = {U{a)-lf. 

The object function for the second dose distribution is given by almost the same formula, but U is 
replaced by the volume function in equation ©. The function U must satisfy several conditions so that 
it is consistent with the dose distributions we chose. First U must be differentiable and decreasing, 
furthermore U(0) > 1 since we assumed that V CT (c? max ) < V a (d c ). A possible choice is U(a) = 10 1 ^ er . 

In the first optimization problem the object function F± is differentiable so there are no difficulties in 
applying the Newton algortihm. However, in the second optimization problem the object function Fi is 
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not twice continuously differentiable at a — 0. Recall that the 1-dimensional Newton iteration for F is 
cr n _|_i = 4>{<J n ) = o~ n — F»(a\ ■ Fixed points o~o of <p correspond to extremal points of F if F" is continuous 
and nonzero at uq . Since the limit from the left of the second derivative of F2 is unbounded, the Newton 
iteration has a fixed point at a — which is not related to an minimum or maximum of F2. But this 
'spurious' fixed point is not stable for the Newton iteration <j). Only if o~ n — the next iterates will also 
be zero. Therefore it does not lead to convergence to a non extremal point of F in general. Although it 
may slow down the convergence of the Newton iterates. For a graphical explanation see figured] 



0~n+l 




Figure 2: In the left figure, the graph of F2 is indicated by a dashed line, F^ by a solid line and F 2 by 
dotted line. The figure on the right shows a magnification of the graph of F' 2 near a = 0. The dashed 
lines indicate how the Newton iterate cr„+i is constructed from a n . 



4 Conclusion 



Our results show that the Radiotherapy treatment planning optimization process based on a quasi- 
Newton method locally converges to an optimum of the object function for most values of the treatment 
parameters. Only when the object function contains a dose- volume constraint (h, v), convergence can not 
be guaranteed. The reason is that the object function is not sufficiently differentiable for values of the 
treatment parameters such that h is a critical value of the dose distribution f a . Although for all other 
values of a, the object function is differentiable. The question remains whether this leads to practical 
consequences for optimization in treatment planning. 

Let us take a closer look at the values of the treatment parameters for which the object function is not 
sufficiently differentiable. In parameter space these values of a are determined by the equation h = / CT (a; CT ), 
see section [3J expressing that h is a critical value of the dose distribution. Geometrically this equation 
defines a hyper surface A in parameter space. Roughly speaking A is m — 1-dimensional so it has measure 
zero. This means that given a starting point for optimization, the Newton iteration scheme generates a 
sequence of points that will lie on A with probability zero. The points on A are spurious fixed points of 
the Newton iteration like in the example of the previous section. The reason is that we can generalize this 
example. First all types of extreme points of the dose distribution, namely minima, saddles and maxima 
lead to the same type of non-differentiability of the volume function. Second, also in higher dimensions 
the second derivative of the object function becomes unbounded at points of insufficient differentiability. 
Third, the spurious fixed points of the Newton iteration scheme are unstable. Our final conclusion is that 
in practice the convergence of the optimization process is most likely only slowed down near A. 



5 Proof of differentiability statements 



In order to make a statement about the differentiability of the object function, we need a few definitions. 
Let £ C R m be an open set of m treatment parameters. Let D be a domain in R 3 representing the 
patient and let R C D represent a tumour region or a healthy organ. For each a 6 S, the dose 
distribution is represented by a function f a : D — > M (here we do not consider f a as a distribution in 
the sense of generalized functions). A function is called smooth if it is infinitely many times continuously 
differentiable. We now assume that / CT is a smooth function of position x £ D and moreover smoothly 
dependent on the parameters a G S. Critical points of / CT are points in D where the gradient of f a 
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vanishes, other points are called regular points. The value f a takes at a critical point is called a critical 
value. If the level set {x £ R \ fo-(x) = h} contains no critical points, h is called a regular value. We now 
state our results on the differentiability of object functions containing EUD constraints and dose-volume 
constraints respectively. 

Theorem 1 Let the object function F : E — > R be given by F(a) = G(E a (a), h), where G : R x R — > R 
is a smooth function and E a (a) is the equivalent uniform dose. Then F is smooth with respect to a . 

Theorem 2 Let the object function F : E — > R be given by F(g) = G(V a (h),v), where G : R x R — * R is 
a smooth function and V a is a volume function for dose distribution f a . Then F is smooth with respect 
to a provided that h is not a critical value of f a . If h is a critical value of f aa , then F is not twice 
continuously differentiable at a = gq. 

The proof of theorem [T] is rather straightforward, but the proof of theorem [3] is more involved. Since 
in the setting of theorem [5] the differentiability of the object function F is completely determined by 
the volume function, the proof is about the latter only. The theorem contains two statements which we 
here state and prove separately. The proofs heavily rely on the results in [4], but due to the presence of 
parameters we have to make several adjustments. Let us start with some definitions. 

Definitions 

1. Let f a : R™ — > R be a positive smooth function for all a £ R m and moreover f a depends smoothly 
on g. For each a £ E, f a vanishes at infinity in the following sense, for each e > there is a 
compact K C R n such that for all x £ R™ \ K, f a (x) < e. We assume that critical points of f a are 
non-degenerate, that is if a; is a critical point of / CT then det(Hess f<j{x)) ^ 0. 

2. Also / : R™ x R m — > R is a smooth function with non-degenerate critical points. Although we 
consider / and f a as different functions, their values are taken to be identical: f(x,a) — f a (x) for 
all x£W l and a £ R m . 

3. We denote the level sets of f a by N h ,„ = {x £ R™ | f a (x) = h}. 

4. The level sets of / are denoted by Af h = {(x, a) £ R" x R m | f(x, a) = h}. 

Remark. Note that critical points of / are also critical points of f a , but the opposite is not necessarily 
true. 

Proof of theorem [lj The function f a is a positive and smooth function, smoothly depending on 
parameters a. Then [f a ] a for a > is also smooth in both variables an parameters. Now it immediately 
follows that E a (a) as defined in equation (fj]) is again a smooth function of a. □ 

The following two propositions essentially cover theorem [H 

Proposition 3 Let f a and f be as defined above. Then V a (h) — vol n {x £ R™ | f a (x) > h} depends 
smoothly on g at a = gq if h is a regular value of f ao . 

Proposition 4 Let f a and V a {h) be as in the previous proposition. Then V a (h) is finitely differentiable 
with respect to g at g = gq if h is a critical value of of f ao . 

Let us begin with the nice situation where h is a regular value of f ao , so that V a (h) smoothly depends 
on g. In the proof of proposition [3] we use a property of the level sets Nh, a of f a . 

Lemma 5 Let f a and f be as defined above and furthermore h be a regular value of f aa . Then an 
open neighbourhood E of Go exists, such that Nh,& is diffeomorphic to Nh,a- for all g £ E and the 
diffeomorphism depends smoothly on g. 



Proof of proposition \3[ The volume V a (h) is given by an integral over a region bounded by Nh.a- 
Since, according to lemma [5J the latter depends smoothly on a at a = gq, the integral and thus V a {h) 
depends smoothly on g at g — gq. □ 
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It remains to prove lemma [SJ 

Proof of lemma [5j Let {a = <j } be a shorthand for {{x,a) G R" x R m | a = cr }- Note that 
Nh,<j = {c = o"o} f~l .A/ft. Since /i is a regular value, Nh jO0 contains no critical points of J^. Thus the 
intersection of {a = <jq} and A/ft is transversal. Then a neighbourhood £ of Co exists such that for all 
(Ji G £, the intersection of {<r = cri} and A/ft is transversal. Thus for all ci G £, Nh, ai is a regular level 
set, consequently there are no critical points of / in A/ft^ = {(x,cr) G R™ x R m | x G Nh, a ,<j G £}. 
Suppose (Ji G S and o\ ^ do then an index i exists such that <j\^ ^ ao^. Now the gradient flow of the 
height function <n on A/ft is a regular flow on Mh.t, thus defining a diffeomorphism $ from A/ft iO . to Nh j(Jl . 

' □ 

Let us now turn to the case where ft. is a critical value of f ao . Then there is at least one critical point 
x ff0 in -/Vft j(To , but for the sake of simplicity we assume that x ao is unique. In the proof of proposition 3] 
we will use a similar construction as in [1] without providing all details. 

Proof of proposition |4l Let us assume that x aa is the unique critical point of / CTo on the critical level 
set Nh,<r Q ■ Since by assumption x ao is a non-degenerate critical point we can apply the implicit function 
theorem to conclude that an open neighbourhood £ of ao exists such that the map E — > R™ : a \— ► x a 
and x a is a critical point of f a is smooth. Now we define a new function g a {x) = f a {x a + x), which is 
a smooth function of x and a G £. Moreover g a has a critical point at x = with critical value f a {x a ). 
Then we have V^(/i) = vol„{a; G W 1 \ g a {x) > h}. Following the construction in [2] we split the latter 
region in a part not containing the point and a remaining small part including 0. The volume of the 
first part depends smoothly on a using the same arguments proving proposition [3] Since V a (h) is the 
sum of both volumes, the differentiability of V a (h) is determined by that of the volume of the second 
part. 

Let us consider this second part contained in a small neighbourhood O of 0. We may assume that O is 
small enough to put g a into the standard form of [4]. In order to keep notation simple we denote the 
standard form again by <? CT , but now we can write g a ix) = g{x) + f a (x a ). The actual standard form is 
■i x p+i> with p + q — n. Note that g no longer depends on a, it has a critical point 
at x = and the critical value is 0. This form can be obtained by a smooth change of coordinates, even 
smoothly depending on a. Here again we need the non-degeneracy of critical point of f a . In the following 
we use the results of [4]. As a function of k, the volume function V\ oc (k) = vol n {a; G O G W 1 \ g(x) > k} 
is only finitely differentiable with respect to k at k = 0. 

Now the volume of the second part is proportional to V\ oc {h — f a (x a )) — vol{a; G O C K™ | g(x) > 
h ~ fa(xa)}- The proportionality factor is of no concern because it originates from smooth coordinate 
changes and depends smoothly on a. Thus we obtain that the differentiability of V a (h) with respect to 
a is determined by that of V\ oc (h — f a {x a )). The latter is only finitely differentiable with respect to a at 
those values of a where h — f G {x a ) = 0. Or, put differently, where h is a critical value of f a . □ 

Remarks 

1. The statements and proofs in this section are valid for any dimension. Therefore we have not 
specialized them to the case n = 3. In this respect finitely differentiable means less than two times 
continuously differentiable. 

2. The critical points of a parameter family of functions will generally be degenerate for certain 
parameter values. However we excluded this possibility because the volume function will be even 
less differentiable at a critical value when the critical point on the level set of this value is degenerate. 
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